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Non-linear transport of electrons in strong electric fields, as typified by dielectric breakdown, is 
re-formulated in terms of the ground-state decay rate originally studied by Schwinger in non-linear 
QED. We discuss the effect of electron interaction on Zener tunneling by comparing the dielectric 
breakdown of the band insulator and the Mott insulator, where the latter is studied by the time- 
dependent density-matrix renormalization group. The relation with the Berry's phase theory of 
polarization is also established. 
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Introduction — While there is a mounting body of in- 
terests in non-linear responses in many-body systems, the 
dielectric breakdown of Mott insulators is conceptually 
interesting in a number of ways. For band insulators the 
dielectric breakdown has been well understood in terms 
of the Zener tunneling[| across the valence and conduc- 
tion bands, which triggers the breakdown. By contrast, 
non-linear properties in strongly correlated electron sys- 
tems should be qualitatively different, since the excita- 
tion gap in the insulating side of Mott's metal-insulator 
transition, arising from the electron-electron repulsion, is 
something totally different from the band gap in its origin 
and nature. Experimentally, Taguchi et al. observed a 
dielectric breakdown by applying strong electric fields to 
typical Mott insulators (quasi-one-dimcnsional cuprates, 
Sr2Cu03, SrCu02)@- The present authors with Arita 
proposed theoretically that the phenomenon may be ex- 
plained in terms of a Zener tunneling for many-body lev- 
els, and obtained the threshold field strength 

Now, there is an important difference between the 
Zener breakdowns in the band and Mott insulators: ex- 
citations (i.e., electrons and holes produced by the field) 
move freely in the former, while they must interact with 
surrounding electrons in the latter and become dissi- 
pated. In other words, while we have only to worry about 
the valence-conduction gap, we have to deal with many 
charge gaps among the many-body levels when the inter- 
action is present. Several authors, including the present 
authors, have shown, with effective models, that a sup- 
pression of quantum tunneling due to the quantum in- 
terference should occur in many-level systems driven by 
external forces (see 0,H,E| and refs. therein). 

The tunneling rate, first obtained by Zener, gives the 
amount of the excitations, and is proportional to the 
leakage current if all the excitations arc absorbed by 
electrodes 0. So the rate is a crucial quantity, but theo- 
retical studies have been quite scant: Even for band insu- 
lators, the quantity is only calculated for systems having 
a simple band dispersion. More importantly, Zener's tun- 
neling rate was based on a one-body WKB approach, so 
we must extend the formalism to many-body systems. 
Specifically, the actual breakdown should be related to 
the scattering between excited states in many-body sys- 
tems as noted above. While this is shown to lead to 



a suppression of the current in a toy (quantum-walk) 
model [(j, we are still badly in need of studies for mi- 
croscopic models, since the existing calculation for the 
Hubbard model|3] was limited to small systems (hence 
to short-time behaviors), while the long-time behavior, 
affected strongly by scattering between excited states, is 
in fact relevant. 

In this Letter, wc propose to determine the tunnel- 
ing rate from the ground- state-to- ground- state transition 
amplitude, whose long-time asymptotic defines the effec- 
tive Lagrangian. The effective-Lagrangian approach was 
evoked by Heisenberg and Euler in their study of quan- 
tum electrodynamics (QED) in strong electric fields^]. 
This was extended by Schwinger to obtain the decay 
rate of the QED vacuumQ. It is well known that 
Schwingcr's electron-positron creation rate can be un- 
derstood by the Landau-Zener mechanism (see e.g. 
What we have done here is the following: (i) We first ex- 
press the effective Lagrangian for band insulators, which 
contains higher-order terms in the Landau-Zener's tun- 
neling probability. Since breakdown of band insulators 
is a well understood subject, the main aim of this sec- 
tion is to introduce the notations for later discussions, 
(ii) Wc then move on to a Mott insulator in strong elec- 
tric fields. The ground-state decay rate is obtained for a 
microscopic (one-dimensional half-filled Hubbard) model 
with the time-dependent density-matrix renormalization 
group [13, 112 • 

From this wc determine the threshold elec- 
tric field to construct the "dielectric breakdown phase 
diagram" . (iii) Finally we comment on an intriguing 
link between Heiscnbcrg-Eulcr-Schwingcr's effective La- 
grangian approach and a recent, Berry's phase approach 
to polarization proposed in Rcfs. |l^llall4 ll5lllq. 

Dielectric Breakdown of a Band Insulator — We start 
with the dielectric breakdown of band insulators in an 
electric field F within the effective-mass approximation. 
For simplicity we take a pair of hyperbolic bands s±(k) = 
iy/V 2 + v 2 k 2 (considered here in d spatial dimensions), 
where V is the band gap, — (+) represents the valence 
(conduction) band, and v the asymptotic slope of the 
dispersion. 

We first obtain the ground-state-to-ground-state tran- 
sition amplitude with the time-dependent gauge in the 
periodic boundary condition. There, a time-dependent 
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AB-flux measured by the flux quantum, 4>{t) = FLr/h 
(with the electronic charge put equal to e = 1, L the 
system size) is introduced to induce an electric field 
F, which makes the Hamiltonian time dependent as 
H{<t>{r)) = Ek,a=±e«( fe + ^(r)e||) C t(fc) Ca (fc). Here 
e || is the unit vector parallel to F, and c) a {k) the cre- 
ation operator with spin indices dropped. If we denote 
the ground state of H((f>) as |0; cj>) and its energy as Eo(cj)), 
the ground-state-to-ground-state transition amplitude is 
defined as 

S(t) = (O;0(r)|fe-^^ ff(< ' i(s))ds |O;0(O))e*^ iJ «(^)) ds ,(l) 

where T stands for the time ordering. The effective La- 
grangian C(F) for the quantum dynamics is defined from 

the asymptotic behavior, E(r) ~ e^ rL C ( F \ The imag- 
inary part of the effective Lagrangian gives the de cay 
rate T(F)/L d = 21m C(F)/h a la Callan-Colemanpjl 
(see also |18|). T(F)/L d gives the rate of the exponen- 
tial decay for an unstable vacuum (ground state), which, 
in the case of Zcncr tunneling, corresponds to the cre- 
ation rate of electrons and holes. Originally, the creation 
rate for Dirac particles was calculated by Schwinger with 
the proper regularization method||. Below, we present 
a simpler derivation which can be extended to general 
band insulators. 

The dynamics of the one-body model can be solved an- 
alytically, since we can cut the model into slices, each of 
which reduces to Landau-Zener's two band model [lHI^ . 
Namely, if we decompose the k vector as (fcj_, fen), where 
fej^ (feii) is the component perpendicular (parallel) to F, 
each slice for a given k± is a copy of Landau-Zener's 
model with a gap Aband(k) = %\/V 2 + v 2 k\. The 
Landau-Zener transition takes place around the level 
anti-crossing on which fey + ^5(/>(r) moves across the Bril- 



louin zone(BZ) in a time interval St 
can be expressed as a transition, 



h/F. The process 



4(fc) y/l -p(k)e- tx(k) cl(k) + >/p(fc)cl(fc), 
cl(k) -Vp(fe)4(fe) + v/1 -p(fc)e lx(fe) cL(fe). (2) 

Here the tunneling probability for each k is given by the 
Landau-Zener(LZ) formula [l9ll2Cj . 



p(k) = exp 



(A band (fc)/2) s 
vF 



(3) 



while the phase x(fe) = —8(k) + j(k) consists of the 
trivial dynamical phase, hQ(k) = J Q T e+(k+^-<j)(s)e^)ds, 
and the Stokes phase |2(ll21|. 



7(fe) = -Im 



ds- 



- 4 (A band (fc)/2) 2 S 



cot(v Fs) — 



vFs 



The Stokes phase, a non-adiabatic extension of Berry's 
geometric phase |23 , depends not only on the topology of 
the path but also on the field strength F [2j| . In terms of 



the fermion operators the ground state is obtained by fill- 
ing the lower band |0; </>) = f^fc c_(k— 3^0eii)|vac), where 
|vac) is the fermion vacuum satisfying c±(fe)|vac) = 0. If 
we assume that excited charges are absorbed by elec- 
trodes [24| we obtain from eqs.Q, J2J) 
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where the dynamical phase 9 cancels the factor 
ei/o E o(tl>(s))ds j n ^ Performing the k integral in 
eq.|@J leads to the ground-state decay rate per volume 
for a d-dimensional hyperbolic band, 



r(F)/L d 



F 



{2-n) d - x h 

CO 



(d-l)/2 



n (d+l)/2 

n=l 



erf 



F 
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To include the spin degeneracy we multiply this by (2s + 
1) for a spin s case. 

If we compare eqs. (@J, (jSJ with Heisenberg-Euler- 
Schwinger's results of the nonlinear responses of the QED 
vacuum in strong electric fields @, II], two significant dif- 
ferences can be noticed. One is the error function ap- 
pearing in eq.JSJ), which is due to the lattice structure 
(with the k integral restricted to the BZ). In the strong- 
field limit (F -> 00), the F^ 1 ^ 2 factor in eq.® is 
cancelled, and the tunneling rate approaches a universal 
function, T(F)/L d -»■ -f In [l - exp (-ttF^/F)] , in- 
dependent of d. Another, quantitative difference appears 
in the Zencr's threshold voltage^; F t b h and = V 2 /va (a: 
lattice constant) for the dielectric breakdown is many or- 
ders smaller than the threshold for the QED instability 

F QED = ™|£l „ 10 16 V/cm 

Dielectric Breakdown of a ID Mott Insulator — Hav- 
ing clarified the one-body case, we now discuss the ef- 
fect of the strong electron correlation on a non-linear 
transport, i.e., the dielectric breakdown in the onc- 
dimensional Mott insulator. We employ the half-filled 
Hubbard model in static electric fields. In the time- 
independent gauge, the Hamiltonian is H + FX with 



H = ~*Eia \ c l+i a Cia + H.cJ +UJ2j »jT"j4> and X = 
Y^j jfij the position operator^) • When F = the 
ground state is a Mott insulator having a many-body 
energy gap so far as U > |26l l2l| . 

After a strong electric field is switched on at r = 0, 
quantum tunneling begins to take place, first from the 
ground state to the lowest excited states, and then 
to higher ones. The tunneling to the lowest excited 
states, which dominates the short-time behavior, has 
been shown to be understood with the Landau-Zener for- 
mula (see eq.(|HJl) 0, but it is the long-time behavior that 
should be relevant to the actual breakdown. 



3 




0.04 
^0.02 

u 



(b) 



^th (U) 



0.0 



0.2 

F/t 



0.4 



FIG. 1: (a) The temporal evolution of the ground-state sur- 
vival probability |H(r)| 2 after the electric field F is switched 
on at t — in the ID half-filled Hubbard model with 
U/t — 3.5, obtained with the time-dependent DMRG for 
L — 50 with the size of the DMRG Hilbert space m — 150 
and the time step dr = 0.02. The dashed line represents 
- In |E(t)| 2 = T(F)t + c for F/t = 0.17. (b) The decay rate 
versus F in the half-filled Hubbard model. Dashed line is a 
fit to eq.0, where F^ ott {U) is the threshold. 



The ground-state-to-ground-statc transition amplitude 
is, in the time-independent gauge, 



S(r) = 



\t(h+FX) 



\0}ei Tl 



(6) 



where we denote the ground state of H as |0) and its 
energy as Eq. The transition amplitude is calculated here 
numerically by first obtaining |0) in the open boundary 
condition with the density-matrix rcnormalization group 
(DMRG) method, and then solving the time-dependent 
Schrodinger equation for the many-body system with the 
time-dependent DMRG [13, El. The ground-state decay 
rate T(F)/L in strong electric fields F is then obtained 
from the asymptotic behavior of E(r). 

Figure ^a) shows the temporal evolution of the 
ground-state survival probability |S(r)| 2 for a system 
with U/t — 3.5. We notice that, as time increases, the 
slope of — ln|S(r)| 2 decreases after an initial stage. The 
slope is proportional to the decay rate, so its decrease im- 
plies a suppression of the tunneling. We can regard this 
as evidence that charge excitations are initially produced 
due to the Zener tunneling, but that, as the population 
of the excitations grows, scattering among the excited 
states become important. This results in the "pair anni- 
hilation" of carriers, which acts to suppress the tunneling 
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FIG. 2: The dielectric breakdown phase diagram on (U, F) 
for the one-dimensional Hubbard model. The symbols are the 
threshold F^ ott (U) obtained by fitting the decay rate T(F) /L 
to eq.@), while the dashed line is the prediction F = F^(U) 
of the Landau-Zener formula eq.JSJ. 



rate. We have determined r(F) from the long-time be- 
havior with a fitting — In |S(r)| 2 = Y(F)t + const. 

The decay rate per length T(F)/L is plotted in 
FiglHTj). where we have varied the system size (L = 
30,50) to check the convergence. T{F)/L is seen to re- 
main vanishingly small until the field strength exceeds a 
threshold. To characterize the threshold F t h(U) for the 
breakdown we can evoke the form obtained above for the 
one-body system (cq.J^Jl). 



IF 

T(F)/L = -—a(U)\n 



_ n FT^U) 



(7) 



(with a factor of 2 recovered for the spin degeneracy). 
The interest here is whether this holds when we replace 
F t b h and with F™ ott (U). The factor a(U) is a parameter 
representing the suppression of the quantum tunneling. 
The dashed line in Fig|T](b) is the fitting for U/t = 3.5, 
where we can see that the fitting, including the essentially 
singular form in F, is surprisingly good, given a small 
number of fitting parameters. The value of a(U) turns 
out to be smaller than unity (taking between 0.77 to 0.55 
as U/t is increased from 2.5 to 5.0). 

In FigH we plot the U dependence of i^ ott . The 
dashed line is the prediction of the Landau-Zener for- 
mula (see eq.Q) 



[A charge (t/)/2] 



(8) 



which was first applied to the dielectric breakdown of 
the half-filled Hubbard model in ref.Q- For the size of 
the Mott (charge) gap we use the Bethe-ansatz result: 



^charge(C^) ~ XI fl 



sinh(2^/^) ^l2fl 123 With v/t = 2. 

So the overall agreement between the threshold for the 
Hubbard model and eq.@ is again confirmed, but to be 
more precise we note the following. Since the scatter- 
ing among charge excitations suppresses tunneling, the 
threshold for the breakdown in interacting systems is ex- 
pected to be larger than the Landau-Zener prediction 



4 



(eq.ljHl)- It was proposed in [6j that the wave function 
in electric fields just below the threshold is localized (on 
energy axis) around the ground state, which physically 
corresponds to a state where the production and anni- 
hilation of excited states are balanced due to the quan- 
tum interference. The F/t = 0.11 result in Fig^a) is 
typical of such states, where — ln|S(r)| 2 first increases 
but soon saturates. When we further increase the field 
strength, the production exceeds the annihilation, which 
leads to an exponential decay of the ground state with 
rate T(F)/L. The system is now metallic in the sense 
that there are carriers that contribute to transport. 

As for the question of whether the dielectric breakdown 
in the non-linear regime may be regarded as a true transi- 
tion, future studies should be needed. The form (eq.®) 
we used for fitting the tunneling rate assumes that the 
transition is a smooth crossover around F ~ F^ ott (U). 
However, the toy (quantum-walk) model calculation in- 
dicates that the breakdown of many-body systems is a 
localization-delocalization transition^. This suggests 
that the tunneling rate and other physical quantities may 
become singular at the breakdown field, in which case 
the I-V characteristics exhibits a jump or a kink. In our 
numerical calculation, while we do find an indication of 
localization for F/t = 0.11 in Fig^a), we observe no 
singularities in the tunneling rate. So, future calcula- 
tions will reveal the nature of the transition around the 
breakdown. 

Finally, we comment on the relation of the effective 
Lagrangian approach with earlier approache s, espec i ally 
the Berry's phase theory of polarization^^ EHU E3- 
Il6[ . There, the ground-state expectation value of the 
twist operator e~ l ^ x , which shifts the phase of electrons 
on site j by — ^-j |l6l| , plays a crucial role. It was revealed 
that the real part of a quantity 

w = -^ln(0| e - 4 T^| ) (9) 



gives the electric polarization, P c \ = — Row [l4|, while 
its imaginary part gives a criterion for metal-insulator 
transition, i.e., D = 47rlmu; is finite in insulators and di- 
vergent in metals [l^j . The effective action in the present 
work is regarded as a non-adiabatic (finite electric field) 
extension of w. To give a more accurate argument, re- 
call that the effective Lagrangian can be expressed as 

C(F) ~ In ((0|e-^ H + F -^|0)e^ Ti50 ) for d = 1. Let 

us set r = h/LF and consider the small F limit. For in- 
sulators we can replace H with the ground-state energy 
Eq to have 

C(F) ~ wF (10) 

in the linear response regime. Thus, the real part of 
Heisenberg-Euler's expression 7] for the non-linear po- 
larization Phe(-F) = — dC gp^ naturally reduces to the 
Berry's phase formula P e \ in the F — ► limit (cf. cq.QJ). 
Its imaginary part, which is related to the decay rate as 
ImPnE(F) = ~^ dr< gp^ L i reduces to — ^ an d gives the 
criterion for the transition, originally proposed for the 
zero field case. 
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After the completion of this work, we came to notice a 
recent paper by Green and Sondhi[28q. where non- linear 
transport near a supcrfluid/Mott transition is studied 
from the point of view of the Schwinger mechanism. 
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